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In this paper we report our systematic calculations of angular momentum / ground 
state probabilities {P{I)) of boson systems with spin / in the presence of random two- 
body interactions. It is found that the P(0) dominance is usually not true for a system 
with an odd number of bosons, while it is valid for an even number of bosons, which 
indicates that the -P(O) dominance is partly connected to the even number of identical 
particles. It is also noticed that the P{Imax)'s of bosons with spin / do not follow the 1/A^ 
(A^ = Z + 1, referring to the number of independent two-body matrix elements) relation. 
The properties of the P(/)'s obtained in boson systems with spin I are discussed. 
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It was discovered a few years ago [1] that the dominance of the 0"^ ground state (0 g.s.) 
of even fermion systems is obtained by diagonahzing a scalar two-body Hamiltonian which 
is randomly determined. Many works have been done in the context of this discovery [2- 
16]. For example, many efforts were made concerning the origin of this observation and 
generic properties of many-body systems in the presence of random interactions. 

In the case of d-boson systems the g.s. probabihty (denoted as -P(O)) was found to 
be close to 0% periodically when the number of d bosons n — 6k±l, where A; is a positive 
integer [10]. This finding of systematic counter examples of g.s. dominance in d-boson 
systems suggests that the breaking of the g.s. dominance might be related to certain 
features associated with bosonic degrees of freedom and thus could be easily seen also for 
bosons with other I. To see whether this is true or not we report in this paper systematic 
calculations of boson systems with spin / in the presence of random interactions. 

The Hamiltonian that we take is as follows: 



where bj is the creation operator of a boson with spin and bi is the time-reversal of the 
annihilation operator bi. The Gl's are two-body interactions between the bosons, and 
are assumed to follow a two-body random ensemble (TBRE), i.e., the G^s are a set of 
random numbers with a distribution function 



All the results in this paper are obtained by using 1000 sets of the above two-body random 
interactions. 

We first look at the case of four bosons with spin /. These systems are interesting 
and worthwhile to study separately, because their fermionic counterparts, namely, four 
fermions in a single-j shell, were studied extensively [7, 10, 11] since these systems are 
the simplest and non-trivial ones to investigate the g.s. dominance of fermions. Fig. 1 
shows a few important angular momentum I ground state probabilities versus I for systems 
with boson number n—A. One sees that the pattern of / g.s. probabilities (P(/)'s) of four 
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bosons with spin I is very similar to that of four fermions in a single-j shell. For instance, 
P(0)'s vs. / stagger with an interval of 6i = 3; P(0)'s arc dominant over other P(/)'s 
except for only two cases: / = 2 and 8; P(/max)'s decrease with /. A new feature is that 
the P(/)'s are considerably large (small) when / is even (odd), i.e., the P(/)'s exhibit an 
even-odd staggering behavior. We should be aware that / = / is possible in the four boson 
case, while in the case of four fermions there are no I—j states. 

The g.s. dominance in the case of four bosons with spin I can be understood based 
on the following observation: there is only one non-zero eigenvalue (denoted as E^^''^) 
among all eigenvalues of / = states if only one of G^, = — 1 and others are zero. This 
observation was proved recently by constructing the 7 = states using a pair basis [16] . 
This means that there is a very large probabihty for 7 = states to give the lowest 
eigenvalue, compared to other / states for which there are many non-zero eigenvalues. 
According to our empirical rule [10], P(/) = f/j/N, where A/} is the number of times that 
the ground state has angular momentum / when one of the two-body matrix elements is 
— 1 and others are zero, and N is the number of two-body matrix elements. The P(0)'s 
are therefore expected to be larger than all the other P(/)'s (with only two exceptions of 
Z = 2 and 8). 

In Fig. 1 one sees that the increase of P(0)'s "coincides" with that of the number 
of 7 = states of four bosons with spin I. Dq^ is 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, • • • 
for l—O, 1, 2, 3, • • •, etc. Although we still can not prove this correlation, below we give 
an argument based on our empirical rule [10]. From the sum rule of diagonal matrix 
elements [17], one obtains that J^lEq^''^ = —\n{n — 1)Dq^ = —6Dq\ where n is the 
number of bosons in the system and n = A here. This relation means that the sum of 
X (-1) increases with -Do"*' suggesting that A/q increases with Dq^ simultaneously. 
Because P(0) = J\fo/N, a regular increase of Dq^ of four bosons with spin I [18] produces 
a regular staggering of the P(0)'s. It is noted that a similar "coincidence" of the P(0)'s 
staggering with increasing number of 7=0 states was firstly observed in [7, 10] for four 
fermions in a single- j shell, but without an explanation or argument. 

We next go to systems with arbitrary particle number n, in order to investigate general 
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features of the -P(/)'s of boson systems with spin / in the presence of random interactions. 
Here a particularly interesting question should be asked concerning the P(/)'s for an odd 
number of particles: which angular momentum / g.s. dominance will appear, / = or 
/ = /? In systems with an odd number of fermions in a single-j shell, there is no /=0 state. 
One expects I = j g.s. probabilities to be large, because one easily relates large j g.s. 
probabilities to large g.s. probabilities for the neighboring even n, and the calculations 
using different types of random interactions showed [7, 10] that this is indeed the case. 
For systems with odd n, however, one may have 7=0 states and thus it is unknown a 
priori whether the I — or I — I g.s. dominance occurs in these systems. 

Fig. 2 shows the -P(0's and P(0)'s in boson systems with I —A and 6, and n running 
as large as possible. For odd n and / = 1 or n = 3 and any odd I, there is no 7 = state; 
for odd n and / = 3, 5, 7 and 9, the 1 = states do not exist unless n > 15, 9, 7, and 
5, respectively. For fifteen bosons with / = 3, nine bosons with / = 5, seven bosons with 
I = 7 and five bosons with / = 9 or 11, P(0) ~ 0% always according to our calculations 
using 1000 sets of a TBRE Hamiltonian. From Fig. 2 and these odd-/ cases we conclude 
that the P(0)'s are in general much less than the corresponding P(/)'s when n is odd. 
On the other hand, the P(0)'s are mostly larger than the P(0's when n is even. This 
result indicates that the g.s. dominance is robust for systems with even n, but not true 
generally if n is odd. In the latter case, we observe that the g.s. dominance is easily lost, 
with a few exceptions. We therefore expect that the g.s. dominance is partly connected 
to the even number of particles. 

It was shown in previous calculations that the g.s. dominance may appear in boson 
systems with odd n, such as sp [8, 9], sd [2, 9, 10], and sdg [10] boson systems. However, 
to a large extent, the very large P(0) therein are associated with the s boson condensation, 
which contributes around 40% of P(0) [9]. In other words, sp and sd systems are very 
special systems in which the s boson condensation produces g.s. dominance when n 
is odd. Without s bosons, the P(0)'s of those systems with odd n would be drastically 
smaller than those with s bosons, and g.s. probabihties for other 7's, such as 7 = Z or 
I — I max, would be much larger. It would be very interesting to carry out systematic 
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calculations of -P(/)'s for both even and odd numbers of bosons with many Vs. 

Now we come to the Imax g-S- probabilities of bosons with spin / interacting by a TBRE 
Hamiltonian. It was shown in [10] that the P(/maa;)'s are associated with the number of 
independent two-body matrix elements, and not sensitive to the particle numbers of the 
system. For fermions in a single-j shell, P{Imax) l/N, where iV = j + |. This relation 
is also applicable to sd and sdg boson systems. Thus it is interesting to check whether 
the PilmaxTs of more general boson systems vs. I show a similar pattern or not. 

Fig. 3 shows the PilmaxV^ vs. I with n ranging from three to six. One sees large 
deviations of the P^maxYs calculated by diagonalizing a TBRE Hamiltonian from those 
predicted by the l/A'" = l/(/-|-l) relation. The P(/max)'s are systematically larger than 
1/A^ (when spin I is small the agreement is rather good, however), and increase with n. 

An argument why the behavior of the P{Imax)^s is different for bosons and fermions 
is as follows. It is known from previous studies [10] that the large P{Imax) comes from 
a gap produced by the pairing interaction Gimax bosons with spin /, and the pairing 
interaction Gj^^^ for fermions in a single-j shell, where Lmax — 2/ and Jmax = 2j — 1, 
respectively. For both cases one may evaluate the gap associated with Gj^^^ or Gimax 
using analytical formulas of {Ej^^^_2 — Ej^^^), because here the state with I^ax {Imax ~ 
2) is found to be the ground (first excited) state if or Gl^^^ is —1. We obtain 

(Ei^a.-2 - Ei^aJ ^ follows: 

2/n- 1 



boson systems 



4/ - 1 ' 



3 105 135 

fermion systems, n = 4 : - + —^rzj-r-. — 37 + 



8 108(4j - 7) 64(4j - 5) 
63 

^128(4j-3)' 

35 2205 5145 

fermion systems, n — o : —— + 



128 2048 (4j - 9) 2048 (4j - 7) 
1785 189 



2048(4i- 5) 2048(4j-3)' 

27 10395 2835 

fermion systems, n = 6 : —— + Tr——7—. — —r + 



128 8192(4j - 11) 1024(4j - 9) 
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4725 45 297 

~4096(4j-7) ~ 256(4j - 5) ~ 8192 (4j -3)' 

(3) 

One easily sees that the gap for bosons with spin I increases regularly with n at an 
interval^ 1/2 if / is large, while that for fermions in a single-j shell is much smaller 
(almost one order) in magnitude and comparable for different n and j. For instance, the 
gap is 0.47, 0.39, 0.35 for n =4, 5 and 6 fermions in a j = 15/2 shell, respectively, while 
the gap is 2.03, 2.56, 3.07 for n = 4, 5, 6 bosons with spin I — 7, respectively. According 
to the empirical rule [10], a relatively larger gap makes the corresponding P{I) larger: 
the larger the gap is, the larger the corresponding P{I) is. Therefore, the PilmaxYs of 
bosons with spin I can be expected to be larger than 1/N = l/(/ + l), which works well 
for fermions and some boson systems with small /'s studied in [10], such as d, sd and sdg 
bosons. 

To summarize, in this paper we reported systematic calculations of angular momentum 
/ ground state probabilities for systems of bosons with spin / in the presence of random 
interactions. First, we found that the g.s. probabilities P(0)'s of four bosons with spin 
I exhibit a pattern similar to that of four fermions in a single-j shell. We explained the 
g.s. dominance and the staggering of the P(0)'s based on a special property of / = 
states when only one of the Gl's are —1 and the other two-body matrix elements are zero. 

Second, we found that the g.s. dominance is usually not true for systems with an 
odd number of bosons with spin / (though the / = states may exist), which indicates 
that the g.s. dominance is partly connected to the even number of particles. 

Third, we found that the Imax g-s. probabilities of bosons with spin / are systematically 
larger than those predicted by the 1/N relation (A^ = 1+1 for bosons with spin I), and that 
P{Imax) increases with boson number n. We presented an argument for this difference 
between bosons and fermions in terms of a large gap for bosons systems with G21 being 
— 1 and others zero. 

One of the authors (ZHAO) would like to thank the Japan Society for the Promotion 
of Science (ID: P01021) for supporting the present work. 
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Fig. 1 The I g.s. probabilities vs. / of four bosons with spin /. The results are 
obtained by 1000 runs of a TBRE Hamiltonian, except that we present the P(0)'s (open 
squares) predicted by a prescription given in Ref. [10], showing that the approach of 
Ref.[10] is reasonably applicable to bosons with spin I. 

Fig. 2 The / = and I = I g.s. probabilities versus n of / = 4 and 6. One sees 
that -P(O) is in general smaller than the corresponding P{1) when n is an odd number, 
indicating that the g.s. dominance is partly associated with the even number of particles. 

Fig. 3. The P(/max)'s vs. / with n ranging from three to six. It is seen that 
the P^ImaxY^ are systematically larger than l/(/ + 1), and increase with n but seem to 
saturate around n — Q. 
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